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Abstract 

Grazing incidence fast atom diffraction (GIFAD, or FAD) has developed as a surface sensitive technique. Compared 
with thermal energies helium diffraction (TEAS or HAS), GIFAD is less sensitive to thermal decoherence but also more 
demanding in terms of surface coherence, the mean distance between defects. Such high quality surfaces can be obtained 
from freshly cleaved crystals or in a molecular beam epitaxy (MBE) chamber where a GIEAD setup has been installed 
allowing in situ operation. Based on recent publications by Atkinson et a/[l[ and Debiossac at a/0, the paper describes 
in detail the basic steps needed to measure the relative intensities of the diffraction spots. Care is taken to outline the 
underlying physical assumptions. 
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1. Introduction 

MBE is the reference technique to produce very high 
quality crystalline surfaces while GIEAD is sensitive to sur¬ 
face coherence over length scales in the range of 1000A 

i i- Recently, combining both techniques has in¬ 
creased the surface sensitivity of MBE diagnostics while 
providing surfaces with exceptional quality for in situ fast 
atom diffraction. GIEAD can be operated easily both 
in growth condition where almost video rate was demon¬ 
strated [l|, and in high resolution mode, before or after 
growth, where its sensitivity poses a challenge to the the¬ 
oretical description [2[. The present paper focuses on the 
data analysis and the associated physics. To support the 
discussion the data presented here are taken from ref.jH 
and [l| correspond to a GaAs surface grown in situ by 
homo-epitaxy and held at a temperature of 570 °G under 
an As 4 overpressure during the measurements. 



Figure 1: Schematic view of a GIFAD arrangement 0, the primary 
beam of helium atoms does not really see individual atoms but are 
reflected by the periodic electronic density of well-aligned rows of 
atoms. 


2. GIFAD and MBE 

GIEAD has been described in several places (see e.g. 
0) and only a brief sketch is given here in fig III and [2j 
Just as HAS, helium atoms are weakly attracted by van 
der Waals forces and strongly repelled by the surface elec¬ 
tronic density. In terms of interaction geometry, as seen on 
figl2l GIEAD is similar to reflection high energy electron 
diffraction (RHEED). In a way GIEAD is to HAS what 
RHEED is to low energy electron diffraction (FEED) as 
illustrated in Table [H 


geometry\ proj. 

electrons 

He atoms 

normal incidence 
grazing incidence 

FEED 

RHEED 

HAS,TEAS 
GIFAD 


Table 1: names of electron and atomic diffraction techniques 

An important practical difference is that keV atoms 
are used which can be detected with high efficiency and 
that the diffraction cone is kinetically compressed allow¬ 
ing the full pattern to fit onto a position sensitive detec¬ 
tor. As a simplication, GIEAD can be seen as a projected 
technique where only the movement perpendicular to the 
surface plane is important 0, 0 • If the helium beam 
is well-aligned with this low index direction forming only 
an angle 0 with the surface plane, the energy of the 
movement normal to the surface is given by = E^sin^O 
with Eq the energy of the primary beam. This energy E^ 
can be tuned between few meV up to few eV with only 
a degree variation as illustrated in fig. [3l Note that the 
wavelength associated with this slow motion is in the 
A range. 
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Figure 2: Schematic view of the MBE chamber with effusion 
cells evaporating gallium and arsenic onto the GaAs(OOl) wafer. A 
He+ ions beam is extracted from a commercial ion source and is 
neutralized before entering the vessel. The atoms scattered by the 
surface are imaged onto a position sensitive detector. 

3. Data analysis 

3.1. primary beam, Laue circle and incidence plane 

We will not consider the nature of the position sensi¬ 
tive detector used to record the diffraction pattern as a 2D 
image. We assume here that the detector is far from the 
surface and perpendicular both to the surface plane and to 
the plane of incidence so that the number of counts in each 
pixel (x^y) corresponds to an intensity map in momentum 
space I{kx,ky). In the present case, one CCD pixel corre¬ 
sponds to a scattering angle of 0.004 deg. or 7 10“^ rad. 
As with standard crystallography two type of information 
can be extracted, the surface lattice unit is reflected in the 
peak spacing while the relative peak intensities depend on 
the scattering amplitudes determined here by the shape of 
the electron density at the surface. The treatment starts 
by a precise determination of the primary beam parame¬ 
ters; its location x^, y^ and its width which is measured 
here by its fwhm = 3.6 pixel or 0.25 mrad, symmetric 
along the x and y directions . This can be achieved before 
or after the diffraction by removing the target surface or, 
during diffraction by leaving a small part of the beam fly¬ 
ing over the surface without interaction, as can be seen as a 
tiny spot at the bottom of figsl3]a),b) and c). In favorable 
cases such as those depicted in flglSl the Laue circle, de¬ 
fined by energy conservation (| kout \ = \kin\) is clearly vis¬ 
ible by interpolation between the elastic diffraction spots, 
and thus its center coordinates Yl and radius are 
easy to pinpoint. However, this does not specify the scat¬ 
tering plane, defined on the detector W a line linking the 
primary beam to the specular beam[2|, [9[. If the direct 
beam is perfectly aligned with the low index crystal axis 
of the surface, the specular spot is easy to identify as the 
symmetry center but the alignment step can be tricky if 
no precautions have been taken 0, 0 • Here we assume 
that the scattering plane is identified as the vertical axis so 
that the horizontal axis of the image is parallel to the sur¬ 
face plane. The goal is now to extract the intensity along 
the Laue circle and to assign, as precisely as possible, the 


intensity of each diffraction spot. 
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Figure 3: Four diffraction patterns recorded for a Eo=400 eV He 
primary beam aligned along the [1-10] direction of the GaAs(OOl) 
surface held at 570 °G corresponding to the /32(2x4) reconstruction 
[2| • The diffractions spots are located on the Laue circle which radius 
is equal to the angle of incidence 6. From a) to d) 6 is increased from 
0.25 deg. to 1.06 deg. allowing a factor close to 20 in normal energy 
Ej_ = Eosin‘^6 between the first and fourth image. 


3.2. Background subtraction 

In figl3]and figlSi) some intensity is present below and 
above the Laue circle. Its contribution increases with the 
angle of incidence, i.e. with the normal energy E± from 
fig St) to figl3]i). This intensity, originating both from 
inelastic scattering 0, and from surface imperfections 
0 [ 0 , is very interesting in itself providing informa¬ 
tion on the surface Debye temperature and on specific 
phonon properties. However, our concern here is only to 
derive the (periodic) topological properties of the surface 
by extracting the elastic diffraction signal, i.e. the inten¬ 
sity associated with the sharp spots sitting exactly on the 
Laue circle. As can be seen in figlU the background to 
be subtracted is far from uniform showing both diffrac¬ 
tion features in the horizontal direction and pronounced 
patches along the vertical direction as well as a quasi- 
nodal structure along oblique directions [l3[ . These struc¬ 
tures originate from quasi-elastic scattering processes 0 
which is closely linked to the elastic scattering 0 . Sim¬ 
ilar patches are indeed present in diffraction charts when 
the diffracted intensities on the Laue circles 0 Q, [lH are 
plotted as a function of the angle of incidence (see also 
figIH]). The reason is that the diffracted intensity can be 
written as an expansion in terms of the number of ex¬ 
changed phonons. Close to the Laue circle, this num- 
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Figure 4: Same data as displayed in figlS]:) and[3]i) but filtered by a 
ID mexican hat function (eq^and insert) in the vertical direction. 
The image looks empty but some of the diffraction spot are already 
saturated is gray scale. 


ber should be limited and the inelastic signal should re¬ 
semble the elastic one in terms of intensity ratio of the 
various diffraction orders [l6|. The fact that this back¬ 
ground looks similar to the intensity on the Lane circle is 
a favorable condition for reasonable corrections. One way 
to estimate this contribution is to consider it as varying 
slowly in the vertical direction so that it can be interpo¬ 
lated from its values above and below the Laue circle by a 
linear interpolation. In practice, the background intensity 
hack{x,y) at any pixel (x^y) on the Laue circle is esti¬ 
mated as hack{x^y) = [I{x^y u) I{x^y — u)]/2 where 
the distance u will be chosen as small as possible but not 
smaller than the experimental resolution. This is equiva¬ 
lent to a double differentiation. To make this subtraction 
easy and not too noisy we use the smoothest possible fil¬ 
ter; a ID doubly differential ’’mexican hat” filter made of 
two gaussian functions (eqUI) (see insert in figH]) allowing 
a straightforward adjustment of the characteristic length 
u (eqU]). 

Fm{y) = 2 * (1) 

Here, a filter with a standard deviation a = 1.4 pixel in 
eqUlis chosen giving a distance u of 3.3 pixels (0.23 mrad) 
between the center of the filter i.e. the positive pole and 
each of the negative poles. This value being close enough 
to the overall resolution taken here as the fwhm of the 
direct beam ; < 75 = 3.6 channels = 0.25 mrad which is still 
significantly smaller than the fastest oscillation rate Vniax • 
This frequency describes how fast a given diffraction order 
switches from bright to dark and then bright again during 
a rocking curve i.e. as the angle of incidence is varied (see 
again figIH]). For a beam energy of 400 eV on the /32(2x4) 
reconstruction of the GaAsfOOli [fR fisj surface Urnax was 
measured at 0.9 0 corresponding here to approxi¬ 

mately 1 mrad or a period P=15 CCD pixels on the detec¬ 
tor. It means that hardly more than four significant data 
points (P/cr 15/3.6) are available to describe this oscilla¬ 
tion. It is therefore very important that the differential 
filter does not degrade the spatial resolution. Figure 
shows that the application of this filter uniformly cancels 
the diffuse background isolating the intensity on the Laue 
circle. The effect of the filter is detailed further in fig. [71 


3.3. Polar-like transformations 

We now need to convert the intensity on the Laue cir¬ 
cle onto a line and this is exactly what a polar transform 
does. However, a standard polar transform centered on 
the Laue circle and associating (kx^ky) (a, \kout\) with 
a = 8iTcsmkx/ky and \kout\ = (^^ + will not pre¬ 

serve the Bragg comb structure kx = n.G. Here n is the 
diffraction order, and G = 27r/a = 0.39A“^ = 6.4 CCD 
pixels, the reciprocal lattice vector associated with the re¬ 
constructed lattice parameter a = 16A. The obvious so¬ 
lution is to keep the kx coordinate intact and to consider 
the transformation (kx^ky) {kx^ |^out|)- To avoid a dual 
assignment, only the particles scattered above the surface 
plane i.e. only the values of ky above the center of the Laue 
circle ky > 0 are considered. For a well defined Laue circle 
as those reported in figs [3] and (H this is enough to gener¬ 
ate a ID intensity distribution shown in fig IT) It should be 
noted that the Laue circle is not alway as easy to identify 
as in figISl this can be due to a limited coherence length 
of the surface defined here as the mean distance between 
crystallographic defects such as, ad-atoms, terraces, or to 
the excitation of surface phonons such as described by the 
Debye-Waller factor adapted to grazing incidence 10. flij. 
This paragraph describes a more robust transformation 
ideally adapted to such common situations. It consists 
in using the primary beam location as a coordinate refer¬ 
ence (0,0). In order to uniquely associate a circle to each 
pixel (kx^ky) another information is needed. We use the 
specular scattering plane to generate a third symmetric 
point {—kx^ky). The specular scattering plane is perpen¬ 
dicular to the surface plane and projects on the detector 
as a line linking the direct beam and the specular beam 
so that only one angle is needed. In practice this angle 
is the global image rotation angle that turns all diffrac¬ 
tion images more or less symmetric with respect to the 
vertical axis. It has to be carefully recalculated for each 
new target surface or each new positioning of the cam¬ 
era but is not supposed to change otherwise. The circle 
associated with each point is now defined as the one en¬ 
compassing itself {kx^ky) , its mirror {—kx^ky) and the 
direct beam (0,0). The radius of this circle defines the ef¬ 
fective momentum keff = ky/2-\-k‘^/ky. The associated 
polar transform writes (kx^ky) {kx^2.keff) as sketched 
in figl5)3) and illustrated in figl6l Note that the factor 2 is 
introduced so that the vertical line is invariant. For points 
located on the Laue circle, both approaches are strictly 
equivalent. Away from the Laue circle, the two trans¬ 
forms differ; when the Laue circle is taken as a reference 
(see figlStt), the surface plane is implicitly taken as a per¬ 
fect reference (all concentric circles stay in this plane) the 
angle of incidence is constant and perfectly defined {\kin\) 
but, away from the Laue circle outgoing scattering angle 
is different (|^owt| 7 ^ |fen|) which poses a serious problem 
if a normal energy and normal wavelength is to be defined 
for diffraction. At variance, in the second approach, the 
direct beam and scattered particle are located on the same 
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circle so that an effective wave vector kef / or angle of in¬ 
cidence can be defined as Oeff = {Oin + Oout)l‘^ he. the 
circle radius in fi^Sb) allowing a quantitative analysis of 
the inelastic diffraction even out of the Lane circle. 



Figure 5: Schematic view of the two options for polar-like transfor¬ 
mations. The reference is indicated by a cross corresponding to the 
origin. In a) the center of the Lane circle is the reference whereas in 
b) the reference is the direct beam. The equation for the respective 
circle radius \kout\ and k^ff are recalled in the insets. 

The underlying assumption is that one has to consider 
the scattering on the surface as a quasi specular scatter¬ 
ing taking place on a portion locally tilted by Oef / — Ogpec 
with respect with the macroscopic plane. This local tilt 
can be understood in classical terms as induced by the 
fact that an assembly of N atoms, each affected by ther¬ 
mal motion usually define a plane which is imperfectly 
aligned. This classical picture considers the surface com¬ 
posed of atoms frozen at their thermally displaced position 
during the scattering. In other words, any given ensemble 
of N atoms give rise to an effective local tilt that can be 
estimated by a linear regression between the given coordi¬ 
nates. 



Figure 6: Same data as in plotted now in (kx,keff) coordi¬ 

nates i.e. after the polar like transformation, a) is without any filter 
showing a complex background structure while b) have been filtered 
in the vertical direction only by the same ID filter as in fig01 

In practice both models, while being strictly identical 
on the Lane circle, are also very similar as long as only 
the region around the Lane circle is concerned. Taking 


the beam as a reference is simply a more pragmatic choice 
with the additional benefit that for each circle an effec¬ 
tive wavelength and wavenumber are associated which can 
be useful for quantitative analysis of inelastic diffraction. 
The result of this transform for figlSt) is depicted in figl6] 
where the vertical direction can now be interpreted as an 
effective wavenumber ke//. The ID vertical mexican hat 
filter described in section [3^ can be applied before or af¬ 
ter the polar transform. Both are almost equivalent in the 
quasi specular region but give rise to a slightly different 
spot shape close to the equatorial plane, i.e. the intercept 
of the surface plane with the detector corresponding to 
the maximum exchange of momentum when the vertical 
momentum kin is entirely transformed into a horizontal 
momentum kout parallel to the surface plane. 

3.4- line profile 



Figure 7: Intensity distribution along the Lane circle in figlSj:) corre¬ 
sponding to the central horizontal lines in figlHi) and b) respectively. 
When the filter is applied,the intensity is reduced but the contrast 
is clearly increased, note that the line shape is affected. 

The intensity distribution on the Lane circle is dis¬ 
played in figHwith and without application of the mex¬ 
ican hat filter. As expected, the intensity is reduced but 
not uniformly. The interesting aspect is that the contrast 
is clearly improved ; some of the lines are almost extin¬ 
guished while no negative intensity is generated (marginally 
at high k). A closer look also shows that the application 
of this vertical filter also affects the horizontal profile of 
the peaks which now have the same gaussian profile as 
the primary beam with exactly the same fwhm. Away 
from the Lane circle the peak profile is clearly Lorentzian. 
It should be noted that without subtraction of the back¬ 
ground, the intensity profile on the Lane circle is poorly 
fitted by gaussians and much better by Lorentzians but 
the present analysis suggests that this is simply due to the 
presence of the background which generates significant in¬ 
tensity between the Bragg peaks. 
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3.5. Diffraction chart 

The relative intensities of all the Bragg peaks contain 
all the information on the scattering of the projectile by 
the surface for a given incidence angle. This is enough 
to compare with detailed calculations but the sensitivity 
is so high that it is not easy to learn from the possible 
mismatch. A diffraction chart i.e. a plot summarizing 
all the intensity profiles associated with the angle of inci¬ 
dence have proven to be much more instructive 0. This 
diffraction chart is simply a rocking curve where all diffrac¬ 
tion orders are plotted. Figure [5] displays such a chart 
in which ve^ distinct patterns can be identified. It has 
been shown [2[ that this pattern can be explained qualita¬ 
tively with only six straight line trajectories illuminating 
the top of the bumps or the bottom of the valleys of the 
potential energy surface. These produce specific path dif¬ 
ferences and the pattern becomes visible just by varying 
the wavelengths, i.e. the phase difference associated with 
these path differences. 



Diffraction order 

Figure 8: Diffraction chart taken from refQ constructed by juxta¬ 
position of successive intensity distributions such as the one in figEI 
at different angle of incidence. A specific pattern is repeated several 
time and depicted in white on the left. 

This simple model [H is summarized here with a brief 
indication of the underlying assumptions. 

-Firstly, the 3D potential energy surface V{x,y,z) can 
be replaced by its 2D average V{y, z) = {V(x, y, z))x where 
X is the low index direction probed (here [110])[20|,[6|. Note 
that it can also be defined as the time average in the pro¬ 
jectile frame V{y,z) = {V{x,y,z))t 

-Secondly, the well documented hard corrugated wall 
approximation considers that most of the momentum trans¬ 
fer occurs very close to the equipotential ID function Z{y) 
defined by V{y,Z{y)) = E± so that straight line trajec¬ 
tories can be considered. As in standard optics when a 
monochromatic light hits a diffraction grating, the diffrac¬ 
tion can be calculated from a simple path difference be¬ 
tween parallel rays. 

-The third simplification qualitatively considers that, 
in a restricted energy range, the shape of the equipotential 
function Z{y) hardly varies with the perpendicular energy 
so that the diffraction chart can now be analyzed as the 
illumination of a fixed profile Z{y) with a variable wave¬ 
length A^. This is valid when E± ^ Eydw where Eydw 


is the typical Van der Waals energy. Note that this situ¬ 
ation is almost never reached with thermal energy helium 
scattering but easily fulfilled with GIFAD. 

- Finally, in the quasi specular region, (particles bounc¬ 
ing almost vertically) the diffraction is considered to be 
dominated by the horizontal (flat) sections of Z{y) i.e. the 
points yi such that dZfyfijdy = 0. These are the the top 
of the hills and the bottom of the valleys (figl9]) i.e. the 
natural points defining the atomic location. 

The whole diffraction chart (inset in figl9j) can now be 
modeled with a reduced set of 6 coordinates fyj^Zj)^ with 
Zj = Z{yj) inside the lattice cell and a straightforward 
analytical formula displayed in eql2] where the six corre¬ 
sponding, phase shifted, unit amplitudes associated with 
any diffraction order n are added and squared. 

Inikin) = aj{n) = nG*yj+2kin*Zj. ( 2 ) 

Note that at F^^=120 meV the effective wavelength 
is only 0.4A; almost ten times smaller that the corruga¬ 
tion amplitude h=3.5 A, the optical analog is therefore a 
grating with a very deep groove. 



0 5 10 15 

y(A) 


Figure 9: The equipotential line Z{y) at 120 meV calculated by 
DFT refQ is replaced by the six point indicated by a cross below 
their coordinates. T is the lattice unit and h=3.5 A is the full cor¬ 
rugation amplitude, inset: Diffraction chart produced by the ray 
tracing model with the six scattering centers. 


The interest of this approach is that it allows an ana¬ 
lytical construction of a diffraction chart such as the one 
in inset of figlU where the topological ingredients can be 
varied at will until a specific feature is identified. Such a 
simple approach is very handy to prepare more elaborate 
description of the surface, full confidence requires confir¬ 
mation by a theoretical construction of the electronic den¬ 
sity profile, and the use of exact diffraction codes such as 
wave-packet propagation 0, UlLclose-couphng [I, 22|, 


or even semi 


i-classical approaches [23|. 


4. GIFAD Images during growth 

One of the main applications of RHEED is the pre¬ 
cise counting of the number of deposited layers. The two 
transformations presented above are simple enough to be 
applied with real time images, but we have focused more 
on the analysis after operation. The analysis of the GIEAD 
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oscillations have been investigated in detail by Atkinson et 
a/[l[ with the homo-epitaxy of GaAs as a reference system. 
In this paper, it was shown that GIFAD oscillations are 
simple, robust and rich. Simple because the phase of the 
observed oscillations is the same for all diffraction orders, 
all primary energies, and angles of incidence. They indi¬ 
cate only the evolution of the surface reflectivity following 
the progression of a new layer. A maximum in a GIFAD 
oscillation indicates the completion of a new layer. The 
oscillations are robust for the same reason that, in layer 
by layer growth, the surface reflectivity oscillates even if 
diffraction features are not resolved. Finally, the GIFAD 
oscillations are rich because the helium scattering com¬ 
prises three different components that were predicted 0 
and observed. These are the elastic diffraction, the in¬ 
elastic or partly coherent diffraction, and the incoherent 
scattering. 



Figure 10: Time evolution of the intensity on the Lane circle during 
growth. It is referred to the intensity before growth. Full circles 
are for the elastic diffraction intensity while open symbol are for 
the total intensity including contributions from elastic, inelastic and 
incoherent scattering. 

The simple mexican hat filter described above was used 
to extract the elastic diffraction intensity and its evolution 
of during growth. Fig|T0] shows that the intensity of the 
elastic component drops to zero at the beginning of a new 
layer whereas the intensity on the Lane circle does not. 
Before growth, i.e. before opening the shutter in front of 
the Ga source the ratio of elastic diffraction, as measured 
with help of the mexican hat filter, is around 36%. This is 
orders of magnitude larger than what is expected from the 
standard Debye-Waller factor for helium scattering from 
a GaAs surface at 570° G. This reduced effective Debye- 
Waller is due to the multiple scattering centers involved in 
grazing angle scattering [l9|,ll0|. The exact behavior of the 
three observed scattering (elastic, quasi-elastic and inco¬ 
herent) intensities is still being investigated. We propose 
here a tentative interpretation in terms of the sensitivity 
to defects. Elastic diffraction probably requires that no 
defect is encountered during close contact with the surface 
while inelastic diffraction could survive as long as the mo¬ 
mentum exchanged with defects is less than a reciprocal 


lattice vector, and incoherent scattering would pose no re¬ 
striction but a is associated with a reduced probability to 
appear exactly on the Lane circle. More work is needed to 
better understand the exact length-scales associated with 
these contributions. As a result a model to fit the relative 
intensities should be possible which would provide data 
with less noise than the double differential filter used here. 
At variance, the advantage of the filter is that it is model 
independent. 

5. Conclusion 

In order to isolate the intensity distribution on the 
Lane circle, two different polar-like transforms have been 
presented taking respectively the center of the Lane circle 
or the direct beam as references. Even in situation with 
a contrasted background is present above and below the 
Lane circle, a simple procedure is described to suppress 
the diffuse background on the Lane circle with a reduced 
statistical noise, allowing an improved contrast in the rela¬ 
tive intensities. The diffraction spots have the same profile 
as the primary beam. 
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